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“When the weather changes, nobody believes the laws of physics have changed. 
Similarly, I don't believe that when the stock market goes into terrible gyrations its rules have 
changed.” 


Benoit Mandelbrot 



RESUMO 


Um dos principais fatores de estudo do mercado de capitais é a discussão a respeito da teoria 
de eficiência de mercado, que no caso diverge em relação ao comportamento do preço da 
maioria dos ativos. Este trabalho tem o intuito de analisar o comportamento do principal 
índice de preços do mercado de bitcoins (BPI) durante o período de julho de 2010 a setembro 
de 2014. Inicialmente será testada a hipótese do passeio aleatório para o BPI. Em seguida 
serão verificadas as correlações de longa data nas séries financeiras temporais utilizando 
como instrumento de análise o expoente de Hurst (H), que inicialmente foi usado para 
calcular correlações em fenômenos naturais e posteriormente sua abrangência alcançou a área 
financeira. O estudo avalia o expoente H através de métodos distintos destacando-se a análise 
R/S e a DFA. Para o cálculo do expoente ao longo do tempo, utiliza-se uma janela móvel de 
90 dias deslocando-se de 10 em 10 dias. Já para o cálculo em diferentes escalas verifica-se, 
para cada dia, o valor do expoente H nos últimos 360, 180 e 90 dias respectivamente. Os 
resultados evidenciaram que o índice BPI apresenta memória longa persistente em 
praticamente todo o período analisado. Além disso, a análise em diferentes escalas indica a 
possibilidade de previsão de eventos turbulentos no índice neste mesmo período. Finalmente 
foi possível comprovar a hipótese de mercados fractais para a série histórica de retornos do 
BPI. 


Palavras-chaves: mercados eficientes, mercado de bitcoins, passeio aleatório, fractais, 
análise de Hurst, análise R/S, análise DFA, autossimilaridade, memória longa. 



ABSTRACT 


One of the main points in the capital market study is the discussion of market efficiency 
theory, in which case differs in relation to the real price behavior of most assets. This work 
aims to analyze the behavior of the Bitcoin Price Index (BPI) during the period from July 
2010 to September 2014. Initially the hypothesis of random walk will be tested for the BPI. 
Then, long period correlations will be observed in financial time series using the Hurst 
exponent (H) as an analytical tool, which was initially used to calculate correlations in natural 
phenomena and then its scope has achieved the financial area. The study calculates the H 
exponent by distinct methods highlighting the R/S and DFA analysis. For the valuation of the 
exponent varying in time, it is used a 90 days moving window dislocating from 10 to 10 days. 
The calculation of the H exponent at different scales analyzes, for each day, the value in the 
last 360, 180 and 90 days respectively. The results show that the BPI index presents persistent 
long memory in almost the whole period of study. Furthermore, the analysis at different scales 
indicates the possibility of predicting turbulent events in the index for the same period. 
Finally, it was possible to prove the hypothesis of fractal markets for the series of BPI returns. 


Keywords: efficient market, bitcoins market, random walk, fractals, Hurst analysis, R/S 
analysis, DFA analysis, self-similarity, long memory. 
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1 INTRODUÇÃO 

O mercado financeiro mundial vem passando, nos últimas décadas, por períodos de 
grandes instabilidades. Os atuais métodos de análise de risco de uma estratégia de 
investimentos como, por exemplo, o Buy-Hold ou alguma fórmula de gerenciamento de 
portfolios como o CAPM desconsideram os pontos causadores das grandes instabilidades nas 
amostras ( outliers ). Observando a distribuição dos log-retomos de preços de importantes 
índices amplamente conhecidos tais como IBOVESPA ou DOW JONES seria relevante 
perguntar: estes pontos muito distantes da média deveriam ser de fato desconsiderados? 
Outros métodos de análise poderiam ser utilizados em amostras de grande volatilidade? 

A possibilidade de prever o comportamento dos preços de ativos financeiros é de 
interesse de várias pessoas que atuam no mercado financeiro. Com tais informações eles 
podem se beneficiar em cada uma das suas áreas de atuação. Porém o processo que produz o 
comportamento de tais valores, ou o processo gerador dos dados, é tão complexo que existem 
muitos estudos para tentar desvendá-los. A complexidade se faz em virtude da gama de 
fatores que influencia na formação dos preços de tais ativos. 

Recentes estudos de análises das oscilações de preços de ativos financeiros levam a 
questionar alguns dos mais importantes postulados da moderna teoria de finanças. Segundo 
Cajueiro et al. (2006), o estudo do comportamento dos preços pode ser feito do ponto de vista 
da modelagem de processos. Procura-se primeiramente identificar um conjunto de fatores que 
interagem entre si para a geração do preço. Após isso, o comportamento de cada um dos 
fatores é modelado individualmente e em conjunto. Tendo em vista a grande dificuldade na 
medição de determinados processos, outras formas de estudo analisam diretamente o 
movimento dos preços, mais especificamente os log-retomos. Neste caso são realizados 
vários testes estatísticos para identificar quais processos poderiam se ajustar de forma 


adequada ao movimento observado. 



14 


Ao se tratar do comportamento do retorno de um ativo financeiro, analisando sua 
linearidade bem como sua distribuição de probabilidade, as teorias que visam a sua 
compreensão divergem constantemente. Fama (1970) ressaltou a Hipótese dos Mercados 
Eficientes (HME), utilizada para descrever um mercado em que informações relevantes são 
corretamente apreendidas pelos preços utilizando a premissa de que todos possuíssem as 
mesmas informações e no mesmo tempo. O autor torna o assunto mais conciso ao dizer que as 
alterações das informações causam sim mudança no preço dos ativos, porém as mudanças 
subsequentes flutuariam aleatoriamente. Assim, sob tal ponto de vista, os preços dos ativos 
não seguiriam um padrão exato, tornando impossível sua previsibilidade ou modelagem. 

De outro modo, há vertentes teóricas que contestam as premissas da HME. Neste caso, 
pressupõe-se que os ativos financeiros representam um sistema caótico. A origem de tal teoria 
é o estudo dos fractais. Addison (1997) define um fractal com uma estrutura que pode ser 
decomposta em partes menores permanecendo inalterada sua geometria. Este fato justifica a 
aplicação da teoria fractal ao mercado de capitais, uma vez que é sabido que o comportamento 
do retorno dos ativos financeiros oscila constantemente. Ainda neste sentido, Mandelbrot 
(2004) afirma que, ao se identificar padrões, é possível haver previsibilidade e que as 
oscilações nos preços dos ativos não são independentes umas das outras. Ao comprovar este 
fato, revela que várias séries financeiras possuem o que denomina “memória”. 

As teorias tradicionais de finanças (pensamento ortodoxo que remonta à tese de 
Bachelier de 1900 e fundamenta modelos como o de Black & Scholes, o de Markowitz, o 
CAPM e suas derivações) falham ao capturar as singularidades presentes em séries temporais 
representativas deste sistema (MANDELBROT, 2004). Estas singularidades podem ser 
apontadas como a não-normalidade dos retornos, a existência de caudas pesadas, as respostas 
assimétricas dos retornos diante de retornos negativos e positivos, a existência de 
heterocedasticidade, entre outros. Desta forma os modelos ortodoxos geram resultados muitas 
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vezes inadequados para descrever a realidade deste sistema complexo marcado pela influência 
de diversas forças e agentes heterogêneos. 

Estudos como os de Mantegna e Stanley (1999) e Bouchaud e Potters (2000) 
evidenciam de que o mercado não é tão eficiente quanto se imagina, seja por violar as 
hipóteses de independência ou as premissas da HME. No mesmo sentido, Lo (1991), 
empiricamente, verificou a não existência de aleatoriedade (passeio aleatório) no mercado 
norte americano, para as décadas de 1960 a 1980. 

Em continuidade aos estudos de memória longa, Mandelbrot (2004) aprimora 
evidências observadas por Hurst (1951), que analisou a dependência duradoura (memória 
longa) nos volumes de altas e baixas do rio Nilo. Analogamente, Mandelbrot acredita ser 
possível de se replicar esta metodologia em outras áreas, tais como, a de finanças. 

Utilizando a linha de raciocínio de Hurst, Peters (1991) aplicou sua estatística aos 
fenômenos econômicos e séries temporais financeiras. Como resultado percebeu que estas 
também se comportavam de modo viesado. Além disso, Mandelbrot e Wallis (1969) 
perceberam que a relação descoberta por Hurst exibia a mesma forma que séries que 
descrevem o movimento browniano fracionário (MANDELBROT e NESS, 1968) e, assim, a 
fórmula de Hurst era aplicável aos fenômenos que seguiam esse movimento. 

De fato, tanto a HME quanto a Hipótese dos Mercados Fractais (HMF) adotam suas 
premissas, realizam seus testes e estudos e auferem seus resultados. O que fica claro é que 
ainda não há uma absoluta certeza o que poderia vir a descartar uma das linhas de estudo. 
Sendo assim, justifica-se a continuidade de estudos relativos ao tema. Além disso, por se 
tratar de assunto relativamente recente, a exploração deste campo tende a provocar evidências 
significativas, relevantes e inéditas. 

Desta forma, conforme mencionado por Santos (2013), modelos que consideram esta 
complexidade dos agentes e estruturas dos mercados componentes do sistema econômico, 
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assim como as conexões não lineares entre os acontecimentos que permeiam este ambiente, 
seriam capazes de capturar as influências dos ciclos de oferta e demanda sobre os preços, que 
implicam em situações de dependências e reversões constantes nos preços que distorcem as 
premissas da HME. 

Muitos dos estudos nesta linha de pesquisa tem se concentrado na avaliação do 
comportamento dos ativos financeiros mais tradicionais, como ações negociadas em bolsa, 
índices de ações, taxas de inflação, de câmbio, mercados agrícolas, entre outros. Ressaltam-se 
os trabalhos de Cajueiro et al. (2006), Malavoglia (2009) e Santos (2013). Todavia, estudos 
direcionados ao comportamento do mercado de bitcoins como sistema complexo são ainda 
restritos ou inexistentes. 

As variáveis econômicas são frequentemente governadas por processos dinâmicos e 
não-lineares e que podem gerar relações de dependência de longo prazo e padrões cíclicos 
não-periódicos com mudanças abruptas de tendências (PETERS, 1991). Para o caso dos 
preços dos bitcoins este comportamento é semelhante tendo em vista a sua intensa 
volatilidade e por seguirem movimentos cíclicos não-periódicos observados na série histórica 
do principal índice de preços do mercado de bitcoins (BPI) durante o período de julho de 
2010 a setembro de 2014. 

A presença destas características acaba gerando séries temporais com singularidades 
como a existência de estrutura fractal e a observação de uma ordem de integração fracionária, 
características que podem ser modeladas por técnicas, como a análise de R/S de Hurst (1951) 
ou a análise DFA (PENG et al., 1994). Conforme ressalta Barkoulas, Labys e Onochie (1997), 
a existência de estrutura fractal nos preços dos ativos levanta questões relacionadas aos testes 
estatísticos tradicionais baseados em métodos padronizados, modelagem econométrica teórica 
e linear dos preços dos ativos, assim como a precificação da eficiência dos mercados e da 


racionalidade. 
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Os trabalhos mais recentes sobre a investigação da dinâmica fractal estão ligados 
diretamente à avaliação da existência de memória de longo prazo nos retornos de ações e 
outros ativos, taxas de câmbio, derivativos financeiros e commodities. A detecção de memória 
de longo prazo indica a existência de uma estrutura de dependência temporal entre as 
observações na série, que podem configurar padrões cíclicos não periódicos que caracterizam 
uma forma especial de dinâmica não linear como a dos fractais. Conforme Lo (1991), esta 
característica pode ser problemática na medida em que é inconsistente com os métodos que 
tratam de processos estocásticos contínuos aplicados aos métodos martingales de precificação 
de contratos futuros e de opções, o que ressalta a necessidade de análise cuidadosa para 
identificação desta característica (SANTOS, 2013). 

Seguindo esta linha de estudo, surge, portanto, um questionamento: O mercado de 
bitcoins possui memória de longo prazo, representando um padrão duradouro? Em outras 
palavras, os impactos sofridos neste mercado hoje, repercutem por longos períodos de tempo? 
Além disso, autossimilaridade em diferentes escalas permite prever eventos turbulentos? 

Dessa forma este trabalho tem por objetivo fazer uma análise do efeito de memória 
longa persistente no mercado de bitcoins , mais precisamente para a série histórica dos log- 
retomos do principal índice de preços deste mercado (BPI) no período de julho de 2010 a 
setembro de 2014. 

Inicialmente será testada a Hipótese dos Mercados Eficientes e em seguida a Hipótese 
dos Mercados Fractais para o BPI, no mesmo período. Em específico, pretende-se avaliar a 
presença da memória longa através, principalmente, da estatística R/S, desenvolvida por Hurst 
(1951) e posteriormente aplicada por Mandelbrot (1972). Em seguida será analisado o 
comportamento do expoente de Hurst ao longo do tempo, introduzido por Muniandy, Lim, 
Murugam (2001), nos índices de ações, Dow Jones, Nasdaq, S&P500, FTSE, Nikkei e 


Ibovespa. 
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Ademais, uma vez confirmada a HMF para este mercado, uma análise do expoente de 
Hurst em diferentes escalas irá verificar, devido à autossimilaridade, a possibilidade de 
previsão de eventos turbulentos causadores das grandes distorções na série histórica dos log- 
retomos do mercado de bitcoins. 

Sendo assim, este trabalho está estruturado em: uma fundamentação teórica, na qual é 
feito um pequeno ensaio sobre o mercado de bitcoins e o expoente de Hurst, respectivamente 
nos capítulos 2 e 3; em seguida, o capítulo 4 trata do processo metodológico utilizado neste 
trabalho; o capítulo 5 mostra os resultados da aplicação dos métodos; e o capítulo 6 as 
conclusões do trabalho. Toda esta sequência de estudo está focada nas flutuações na série 
histórica dos log-retornos do mercado de bitcoins. 
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2 O MERCADO DE BITCOINS 

Este capítulo aborda basicamente a questão do funcionamento do mercado de bitcoins 
com algumas de suas características. Vale ressaltar que se trata apenas de um ensaio sobre o 
tema, haja visto que este assunto apresenta uma grande complexidade que não comporta ser 
analisada nesta dissertação. 

Inicialmente será abordada uma visão geral da estrutura bitcoin e em seguida citar 
alguns acontecimentos marcantes que ocorreram ao longo de sua história. Nas seções 
seguintes serão feitas explanações sobre os agentes deste mercado, suas regras e riscos. E para 
finalizar, será descrito como é composto o principal de preços deste mercado, o BPI. 

2.1 Uma breve definição do bitcoin 

Bitcoin geralmente é definido como um tipo de moeda virtual, "cripto-moeda", surgido 
via Internet, supercomputadores e a vontade de indivíduos buscando novas formas de troca 
monetária (BITCOIN, 2013). 

Bitcoins compartilham muitas semelhanças com outras moedas, mas o mais 
importante é que mais e mais comerciantes, varejistas e indivíduos aceitam bitcoins como 
pagamento. É possível comprar alimentos e roupas com bitcoins , subscrever a serviços on¬ 
line, ou mesmo fazer compras em supermercados. 

No entanto, bitcoins também são muito diferentes das moedas tradicionais. Ao 
contrário de reais, dólares ou euros, bitcoins não são apoiados por um governo. Eles são uma 
forma descentralizada de dinheiro e não estão ligados a um sistema de banco central ou 
autoridade emissora. Além disto esta moeda existe em um mundo online conduzido via 
algoritmos matemáticos e protocolos de criptografia inteligentes (DACHIS, 2013). 

É possível usar bitcoins para todos os tipos de transações reais. Para fazer isso, basta 
comprar bitcoins da maneira que quiser, seja através de cartão de crédito, conta bancária ou 
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até mesmo anonimamente com dinheiro. Em seguida, os bitcoins são transferidos diretamente 
em uma conta biícoin, aberta pelo próprio usuário via internet em sites de qualquer bolsa de 
bitcoins , seja no Brasil ou em outro país. Assim pode-se enviar e receber pagamentos 
diretamente a um comprador ou vendedor, sem a necessidade de um típico intermediário, 
como um banco ou empresa de cartão de crédito. 

Ao pular o intermediário na transação, paga-se muito menos em taxas associadas. 
Cada uma das partes no acordo também pode manter um nível muito mais elevado de 
anonimato, que tem os prós e contras de todos os envolvidos. 

Pagar ou receber bitcoins é tão fácil como enviar um e-mail, via computador ou 
smartphone. Essa simplicidade esconde o fato de que há um complexo algoritmo de 
matemática para proteger todas essas transações e manter a legitimidade e segurança. 

2.2 O PASSADO DO BITCOIN 

A gênese do biícoin é uma lenda da Internet. Em novembro de 2008, uma pessoa (ou 
pessoas) que trabalham sob o pseudônimo Satoshi Nakamoto publicaram um documento que 
define a viabilidade do conceito de biícoin. Nakamoto mencionou a crise de 2008, bem como 
as falhas de divisas e corrupção de sistemas bancários existentes apoiados pelo governo, como 
um fator motivador para inventar uma nova moeda (BITCOIN, 2013). 

Em 2009 Nakamoto lançou o primeiro aplicativo biícoin , e também realizou a 
“mineração” dos primeiros bitcoins para circulação. Então era só uma questão de espalhar o 
conceito da nova moeda para o mundo. 

Para usar bitcoins , é necessário de um software de carteira de bitcoins , que criptografa 
e mantém atualizado o saldo de bitcoins em um computador ou smarphone. Para começar, 
deve-se primeiro baixar e instalar o software e, em seguida, preencher a carteira com bitcoins 
usando conta bancária, cartão de crédito ou outra forma de pagamento. Este processo de 
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compra e venda é realizado via as bolsas de bitcoins como, por exemplo, a Mercado Bitcoin 
no caso do Brasil ou a Bitstamp e BTC-e no caso de bolsas estrangeiras. É possível também 
transferir bitcoins entre pessoas via software desde que o processo seja monitorado (ou 
“ rnined ”) pelos responsáveis a fim de evitar fraudes (DACHIS, 2013). 

Uma vez comprados e/ou transferidos é uma questão de encontrar um fornecedor que 
aceite bitcoins como pagamento. Quando os bitcoins foram lançados, o volume trafegado era 
muito pequeno, mas nos dias de hoje há muitos estabelecimentos que aceitam este tipo de 
moeda. Isso inclui restaurantes, lojas de roupas, médicos dentre outros. Algumas pessoas 
ainda utilizam bitcoins para aluguel de imóveis e compra de veículos. 

Vale a pena ressaltar as principais empresas dos Estados Unidos que aceitam bitcoins 
como moeda: Amazon, Netflix, Delta, Starbucks, Whole Foods Market, ebaY, HP dentre 
outras (PETRONZIO, 2013). 

2.3 O MERCADO DE BITCOINS COMO SISTEMA COMPLEXO 

Bitcoins não são como ouro e diamante que no caso são descobertos e comercializados 
de forma tangível. Eles também não são como o dinheiro em cédulas. Nenhuma autoridade 
central imprime cédulas de bitcoins tangíveis de difícil falsificação para circulação. Em vez 
disso, os bitcoins dependem inteiramente uma rede de informática descentralizada. 

Todo o sistema de bitcoins é executado em uma rede P2P (peer-to-peer ). Esta 
arquitetura P2P é semelhante à das redes de compartilhamento de arquivos onde pessoas 
distribuem gratuitamente os dados de todos os tipos, incluindo música com direitos autorais, 
filmes e muito mais. É um sistema flexível análogo aos uTorrents. 

Em outras palavras, não há um servidor central executando todos os processos 
relacionados com o bitcoin. Em vez disso, o computador de cada usuário bitcoin é parte da 
rede, compartilhando coletivamente a carga computacional de gerar bitcoins e efetuar e 
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validar suas transações. É esta natureza descentralizada que faz com que o bitcoin não sofra 
intervenções de instituições governamentais, e seja livre de regulação e fiscalização 
(BITCOIN, 2013). 

Antes de utilizar bitcoins as moedas precisam ser analisadas por um processo chamado 
de mineração bitcoin , ou mining. Qualquer computador pode começar a mineração de bitcoins 
usando um aplicativo gratuito de mineração. O mining requer antes que toda a rede de 
participantes de bitcoins processe um determinado algoritmo ou realize um cálculo numérico 
complexo. A cada dez minutos é lançado um novo desafio e o usuário que o desvendar 
garante esta janela de tempo para realizar o mining e, uma vez concluído trabalho, será 
recompensado com a moeda. Assim são criados os bitcoins (NAKAMOTO, 2008). 

Muitas pessoas realizam grandes investimentos em supercomputadores para desvendar 
os algoritmos e garantir a janela de tempo para o mining. Aquele que estiver atualizado, tanto 
em nível de hardware quanto software, garantirá a janela para o trabalho e, 
consequentemente, terá mais lucro (ANDERSON, 2013). 

A quantidade exata de trabalho necessária para a análise das transações não é fixa. A 
rede se ajusta à carga de trabalho de modo que o número de bitcoins emitidos cresça a uma 
velocidade constante, pré-determinada. O processo de emissão continua até que o número de 
bitcoins em circulação atinja o seu número final, definido em 21 milhões. Atualmente, o 
processo de mineração cria 25 bitcoins a cada 10 minutos. A cada quatro anos, o número de 
moedas que podem sofrer o mining é reduzido para metade, até que o limite coberto de 
moedas seja alcançado no ano de 2140. Após esse momento, o número de bitcoins em 
circulação se tornará estático. Isto provavelmente para se ter um controle da emissão da 
moeda em circulação bem como sua valorização (DACHIS, 2013). 

Conforme mencionado anteriormente, é possível acumular bitcoins em uma carteira 
digital. Ao comprar, enviar ou receber moedas, elas são verificadas por uma assinatura digital, 
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chamada de chave de criptografia pública, o que impede a falsificação das moedas e as 
tornam torna reconhecíveis à rede. De certa forma, a rede compartilha a contabilidade de 
todos os bitcoins (o chamado blockchain ), que fornece uma maneira confiável e redundante 
de manter o número de bitcoins em circulação. O blockchain é basicamente um arquivo de 
tamanho considerável (aproximadamente 20GB) e qualquer um pode baixa-lo a qualquer 
momento (BITCOIN, 2013). 

Tudo isso funciona graças ao código open-source do bitcoin (ou seja, visível para 
todos). O software open-source é comumente usado por programadores que se opõem à 
especulação corporativa e controle. Qualquer programador experiente pode ver como o 
algoritmo do bitcoin funciona, e isso não é uma ameaça. Isto porque não é o código que 
protege transações, mas é o fato da rede ser compartilhada que garante a legitimidade de cada 
transferência, pois é de interesse de todos. Uma transação inválida por ser aceita por um 
minerador inescrupuloso, mas quando o falso bitcoin passar para um próximo estágio, seja 
compra, venda ou saque, o novo minerador descobrirá que o anterior foi desonesto e desta 
forma este será excluído da rede. 

2.4 A IDENTIDADE DA MOEDA 

Quando se guarda bitcoins em uma carteira digital, pode-se simplesmente deixa-los lá 
e esperar que o preço dos mesmos se valorize, ou pode-se trocá-los para em moeda local, seja 
dólar, euro, real dentre outras. Caso os mesmos permaneçam armazenados em um 
computador ou smartphone , é imperativo lembrar que não há nenhuma empresa ou entidade 
com um backup dos bitcoins (DACHIS, 2013). 

Isso significa que é fundamental ter um backup da carteira. É melhor guardar esse 
registro em um dispositivo, como uma unidade de memória flash, ou pendrive, que seja 
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possível guardar em um local seguro. Caso contrário, caso o disco rígido do computador 
corrompa, ou caso o smartphone se perca, todos os bitcoins desta carteira serão perdidos. 

As operações com bitcoins são irreversíveis e geralmente ocorrem em um curto espaço 
de tempo. No entanto, pelo fato do processo de autenticação compartilhar dados referentes à 
transação com toda a rede, às vezes pode-se esperar alguns minutos antes que um pagamento 
seja concluído (BITCOIN, 2013). 

Como não há regulamentação nacional para bitcoins, é possível transferi-los para 
dentro ou fora de qualquer país evitando assim as taxas geralmente cobradas pelos bancos e 
agências de câmbio. Além disto, como nenhuma autoridade governamental fiscaliza as 
transações, não há limites de valores a serem transferidos (MELENDEZ, 2013). 

É neste ponto que surgem várias perguntas sobre a legitimidade de bitcoins, ou seja, a 
sua identidade. Como uma moeda pode surgir da noite para o dia pela internet e possuir um 
valor financeiro? Os economistas podem oferecer uma explicação longa e filosófica sobre a 
história do dinheiro, mas a resposta é simples e curta: toda moeda possui valor porque existe 
uma crença na sociedade a respeito deste valor (KRUGMAN. 2013). 

O bitcoin não é diferente nesse aspecto. A moeda tem sido abraçada por especuladores 
financeiros e pessoas que simplesmente não confiam mais sistemas bancários controlados 
pelo governo. Há uma grande confiança no algoritmo de criptografia do sistema bitcoin, 
sendo esta contagiante, dando ainda mais legitimidade a esta moeda virtual. 

2.5 Um mercado de alto risco 

Conforme mencionado anteriormente, o bitcoin está ganhando mais legitimidade 
desde seu lançamento, em julho de 2010. Mas, há a questão da confiança desta como uma 
moeda. E é neste ponto que esta nova moeda se torna espantosa. A especulação desenfreada 
transformou a valorização dos bitcoins em uma montanha-russa. Em 2010, uma unidade de 
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bitcoin valia apenas alguns centavos, mas aumentou de forma constante para cerca de US$30 
em junho de 2011. Perto do final desse ano, porém, despencou violentamente para apenas 
cerca de US$2. 

Em abril de 2013 teve uma recuperação notável chegando a mais de US$140. Em 
seguida, mais uma vez, foi para abaixo de $100. Mas no final deste ano que atingiu seu topo: 
quase US$1150 em dezembro de 2013. Hoje em dia está na faixa de US$500 (CHAPPELL, 
2013). 

Para muitos observadores, esta volatilidade espantosa despertou uma série de bolhas 
especulativas. Alguns economistas alertam que o sistema bitcoin é uma bolha esperando para 
estourar. Há uma crença de que esta moeda foi construída apenas para os especuladores 
(SIDEL, 2013). 

Há também o fato de que o câmbio on-line continua sendo o motivo da maioria das 
transferências de bitcoins. As trocas são necessárias, de fato, pois ajudam a converter bitcoins 
para moeda local. No entanto, em um sistema projetado para ser descentralizado, essas trocas 
oferecem um alvo saboroso para os órgãos regulatórios do governo e hackers mal- 
intencionados. 

Hackers, por exemplo, podem, ainda, não ser capazes de explorar o algoritmo do 
sistema bitcoin, a fim de criar moedas falsas ou transações fraudulentas. Mas eles podem 
definitivamente atacar as bolsas de bitcoins , que são tão vulneráveis como qualquer outro tipo 
de site. Em um incidente amplamente divulgado, os hackers lançaram um DDoS ("ataque de 
negação de serviço distribuído") na famosa bolsa “Mt. Gox”, tomando o site inacessível para 
praticamente todos os usuários. Repentinamente o chãos se alastrou e valor do bitcoin 
despencou (MCMILLAN, 2013). 

As próprias bolsas são arriscadas de certa forma. Muitos usuários armazenam seus 
bitcoins em suas contas cadastradas nas bolsas sob o pressuposto de que esses sites são o 
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método de armazenamento mais seguro. Mas quase metade das bolsas de bitcoins 
desapareceu, levando consigo as reservas de muitos usuários. Por isto é recomendável que a 
custódia permaneça com o usuário, seja em seu computador ou smartphone (THE 
ECONOMIST, 2013). 

Desta forma o bitcoin é, definitivamente, um ativo de extremo alto risco. E como será 
descrito em breve, há outras desvantagens potenciais para este sistema. 

2.6 Atividades extralegais 

O relativo anonimato do bitcoin poderia, de fato, permitir a evasão fiscal. 
Tradicionalmente, esconder ativos em alguma conta no exterior leva algum trabalho braçal 
real. Com o bitcoin , porém, é possível esconder milhões de dólares em uma carteira digital, 
com segurança fora da vista de todos os cobradores de impostos, e esse processo poderiam ser 
concluído em apenas uma tarde (THE ECONOMIST, 2012). 

Governos e cidadãos conscientes, obviamente, se preocupam com estas questões. As 
autoridades não querem deixar de receber impostos, bem como as agências reguladoras, 
compreensivelmente, não estão satisfeitas com o fato das pessoas encontrarem novas 
maneiras de contornar as leis. 

Apesar deste tipo de questões legais ou éticas, e mesmo que o seu valor seja 
totalmente imprevisível, o bitcoin tem resistido. Isso deixa uma série de economistas e 
consumidores curiosos sobre o futuro da moeda. 

2.7 Bitcoin é apenas o começo 

Os picos e vales no valor bitcoin são sinais de alerta claros de que esta nova “cripto- 
moeda” é tudo menos estável. Até que se comprove que este moeda seja utilizada para a 
compra de bens e serviços, em vez de apenas especulações, é difícil dizer se bitcoins 


continuarão existindo. 
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Grandes investimentos são realizados em supercomputadores dedicados 
exclusivamente à mineração de bitcoins. A competição entre os mineradores é tão intensa que 
em determinados momentos o valor gasto em componentes eletrônicos é maior que o prêmio 
recebido pelo processo de mineração (BROWN, 2013). 

Há também uma série de questões não respondidas sobre a natureza da própria moeda 
virtual. Ainda não se sabe se o mundo está realmente pronto para mudar do papel moeda para 
uma versão puramente online. 

Ainda se tratado de cripto-moedas, o bitcoin não é a única no mercado. Há uma 
abundância de concorrentes no mercado de moedas virtuais. Incluem PPcoin, Litecoin, 
Ripple, Freicoin e Namecoin, entre outros. Cada um desses sistemas alega ter pontos fortes 
que outras moedas não podem igualar (SIMONITE, 2013). 

Tratando-se de bitcoins , dólares, reais ou euros, a economia global exige confiar em 
um sistema de moeda que faça sentido e seja consistente. É possível que, nos próximos anos, 
uma moeda digital possa transmitir confiança suficiente para começar a substituir as moedas 
mais tradicionais. Até então, no entanto, o fenômeno bitcoin pode desvanecer-se, pouco a 
pouco, no éter Internet. 

2.8 O ÍNDICE COINDESK 

Lançado formalmente como um índice em setembro de 2013, o índice de Preços ao 
Bitcoin “CoinDesk” (BPI) representa uma média de preços bitcoin de todas as principais 
bolsas globais que atendam a critérios especificados pelo BPI detalhados a seguir. 

Destina-se a servir de referência preço padrão de varejo para os participantes do setor 
e profissionais de contabilidade. Sua fórmula é flexível o suficiente para manter apenas as 
taxas de câmbio que respondam a certos critérios mínimos para determinação de preço e 


validade (COINDESK, 2013). 
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De acordo com os desenvolvedores do CoinDesk, os critérios estão sujeitos a 
modificações com base na evolução do mercado e condições regulamentares que pretendam 
manter as eventuais revisões a fórmula em uma programação mensal previsível. Isto irá 
aumentar a credibilidade do índice, pois os usuários terão alguma antecedência em relação às 
atualizações de fórmulas e fundamentos associados a essas atualizações. 

A criação de um novo índice pode demandar um esforço complicado porque ele deve 
apelar para o menor denominador comum para ser amplamente aplicável. No entanto, deve 
ser seletivo o suficiente para maximizar a descoberta de preços com base no mercado. 

Desde o dia 01 de julho de 2013, as seguintes bolsas de bitcoin cumprem os critérios 
e, portanto, estão incluídas no cálculo do BPI em dólares: 

• Bitfinex 

• Bitstamp 

• BTC-e 

• LakeBTC 

Os principais critérios especificados pelo BPI são (COINDESK, 2013): 

• O CoinDesk BPI é uma média simples dos preços das bolsas que o compõem. 

• O BPI é atualizado a cada 60 segundos. 

• A nova série histórica do índice começa no dia 01 de julho de 2013. 

• Dados históricos antigos são obtidos via bolsa Mt. Gox (desde 18 de julho de 

2010 ). 

• O timestamp do BPI (abertura, fechamento, máximo, mínimo) é baseado no 
Tempo Universal Coordenado (UTC). 

Com relação à bolsa Mt. Gox já mencionada anteriormente, seu preço, de julho de 
2010 a julho de 2013, foi a referência para formação do índice. Entretanto muitas empresas 


tiveram problemas ao considera-la como benchmark de suas operações visto que em 
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determinados momentos a plataforma Gox não estava disponível, visto que foi vítima de 
ataques de hackers (MCMILLAN, 2013). 

Devido aos constantes atrasos nas transferências bancárias da Mt. Gox , um spread de 
cerca de 10% surgia entre ela e outras bolsas por vários meses, causando distorções nas 
conversões em tempo real do índice. Em determinados momentos chegou a 20%. 
Consequentemente a Mt. Gox deixou de fazer parte BPI. 

A decisão de aplicar uma média simples para o BPI, ao contrário da média ponderada 
por volume, se deve ao fato do mercado bitcoin atualmente ter pouca liquidez regional. 

Como agora o volume de negociação favorece determinadas regiões, uma abordagem 
ponderada por volume não agiria como um indicador global adequado, pois cada bolsa bitcoin 
internacional está igualmente disponível para todos os participantes do comércio nacional 
(COINDESK, 2013). 

A média simples não favorece uma bolsa regional com alto volume e garante que o 
BPI seja significativo para o maior número de participantes do mercado. Além disso, a 
abordagem média simples minimiza o impacto das irregularidades no volume e, 
acidentalmente, excluindo uma bolsa. 

À medida que liquidez geral melhore e o número de bolsas no âmbito global aumente, 
o impacto das variações regionais deve diminuir e uma abordagem ponderada por volume 


poderá ser mais apropriada. 
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3 REVISÃO DE LITERATURA 

Este capítulo aborda basicamente algumas informações sobre o expoente Hurst e 
alguns métodos para determiná-lo. Inicialmente será feita uma breve descrição da busca por 
modelos de preços de ativos financeiros com especial ênfase em leis de potência (invariância 
de escala), discutindo as principais contribuições de Bachelier, Mandelbrot e Peters. Em 
seguida serão feitos pequenos comentários sobre memória de longo prazo. Na sequência, 
serão abordadas informações sobre o expoente Hurst e alguns métodos para determiná-lo 
focando principalmente no método R/S, tanto para todo o período da série quanto ao longo do 
tempo. Finalmente serão feitos pequenos comentários sobre persistência e grau de eficiência 
de mercados extraídos do expoente de Hurst. 

3.1 OS ESTUDOS DE BACHELIER, MANDELBROT E PETERS 

O primeiro modelo formal para a dinâmica de preços foi apresentado por Bachelier 
(1900). Seu processo de variação de preço, em essência, define a matemática do movimento 
browniano antes que Einstein e Wiener replicassem seus resultados em 1905 e 1923 no 
contexto das partículas físicas e, em particular, gerando uma distribuição normal (ou seja 
Gaussiana) onde a variância aumenta proporcionalmente com o tempo. Um pressuposto 
fundamental da abordagem de Bachelier é que sucessivas mudanças de preços são 
independentes. Sua tese permaneceu desconhecida para o mundo financeiro até Osbome 
(1959), que apesar de não fazer nenhuma referência ao trabalho de Bachelier, redescobriu 
movimento browniano como um modelo plausível para mudanças de preços de segurança 
(CLARKSON, 2002). 

Em 1963, Benoit Mandelbrot apontou que as caudas das distribuições de preços são 
muito mais pesadas do que as de distribuições normais. Ele observou também que a suposição 
de independência dos dados dos modelos anteriores não reflete totalmente a realidade, 
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dizendo que mudanças longas eram seguidas de mudanças longas assim como mudanças 
curtas eram seguidas de mudanças curtas. Mandelbrot (1963) deparou com o trabalho de 
Hurst (1951, 1955), que tratou de uma forte dependência em eventos naturais de um ano para 
outro e desenvolveu o expoente Hurst “H” como uma medida estatística robusta de 
dependência. No novo modelo de Mandelbrot, “movimento browniano fracionário”, o qual é 
descrito em detalhe em Mandelbrot & Van Ness (1968), foi definida por uma equação que 
incorpora o expoente Hurst H. Muitos economistas, particularmente Cootner (1964), foram 
críticos do trabalho de Mandelbrot, principalmente porque - se ele estava correto sobre 
distribuições normais sendo incoerentes com a realidade - a maioria dos trabalhos estatísticos 
já publicados e consolidados, especialmente tratando-se do CAPM e da HME, seriam 
inválidos. 

Devido à complexidade matemática do movimento browniano fracionário, e em parte 
por causa da conclusão em Lo (1991) de que distribuições padrão, como a normal por 
exemplo, podem dar uma representação adequada da realidade, os dois modelos sugeridos por 
Mandelbrot não conseguiram fazer um grande impacto na teoria de finanças, e ele 
essencialmente deixou a cena financeira para perseguir outros interesses, tais como geometria 
fractal. No entanto, em seu “The Fractal Geometry of Nature”, Mandelbrot (1982) comentou 
sobre o que considerou como as metodologias estatísticas "suicidas", que eram padrão na 
teoria de finanças: 


Diante de um teste estatístico que rejeita a hipótese 
browniana que as mudanças de preços são Gaussianas, 
o economista pode tentar uma modificação após outra 
até que o teste seja fraudado. Uma correção popular é a 
‘censura’, hipocritamente chamada de 'rejeição de 
outliers'. Distinguem-se as ‘pequenas’ mudanças de 
preços comuns a partir das grandes mudanças... Os 
primeiros são vistos como aleatórios e Gaussianos, e 
tesouros de ingenuidade são dedicados a eles... Os 
últimos são tratados separadamente, como ‘não- 
estocásticos’. 



32 


Em 1999, Mandelbrot produziu um artigo (reportagem de capa da edição da “Scientific 
American ” de Fevereiro de 1999), no qual foram utilizadas analogias para destacar a natureza 
temerária dos modelos de risco padrão que assumiam distribuições normais independentes. 
Foi apontado que uma representação mais realista das flutuações do mercado, ou seja, 
movimento browniano fracionário em escala de tempo multi-fractal, já existia. 

O gestor de investimentos norte-americano, Edgar Peters, palestrante frequente e 
seminarista para o público de profissionais da área de investimento, não só acreditava que os 
ensinamentos da teoria de finanças dominante eram fatalmente falhos, mas também viu um 
grande potencial em novas metodologias que levavam em conta o fato do expoente Hurst ser 
significativamente em excesso ao valor de referência de 0,5 do movimento browniano. 

Em seu artigo de 1989 do “Financial Analysts Journal ”, Peters (1989) estima valores 
do expoente de Hurst para o índice S&P500. Em seguida, estabelece, com uma ênfase maior 
na teoria do caos do que na análise da auto similaridade em escala, suas razões para rejeitar as 
principais hipóteses da tradicional teoria de finanças (PETERS, 1991). Este último trabalho, o 
que é, essencialmente, menos complexo matematicamente em seu corpo, mas com apêndices 
altamente técnicos e pertinentes, conseguiu trazer as obras de Hurst e Mandelbrot para um 
público muito mais vasto do mercado financeiro. Posteriormente Peters (1994) expande seus 
trabalhos anteriores e, em especial, inclui detalhes muito mais práticos e complexos na análise 
fractal e movimento browniano fracionário. 

De acordo com Peters (1991), Mandelbrot não estabeleceu em seus trabalhos uma 
definição formal e precisa sobre uma estrutura fractal, apesar do desenvolvimento de extensas 
pesquisas sobre estas estruturas em diferentes áreas do conhecimento. Desta forma, a 
determinação de uma estrutura fractal implica na identificação de um conjunto de regras que 
podem ser percebidas, como ressalta Peters (1994), por meio da investigação quanto às 


características comuns e mensuráveis a estas, tais como a autos similaridade e a dimensão 
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fractal. A autossimilaridade, em matemática, indica que parte do objeto ou estrutura em 
análise é semelhante e está de alguma forma relacionada com o todo ao qual ela pertence. Isso 
implica que as partes de um sistema, em diferentes escalas espaciais ou temporais, são 
semelhantes ao sistema como um todo (FEDER, 1988). 

3.2 A MEMÓRIA DE LONGO PRAZO 


Harold Edwin Hurst (1880 - 1978) passou quase toda a sua carreira profissional como 
um hidrólogo no Egito lutando com o problema do controle de reservatórios. Como Peters 
(1991) observa: 


Um reservatório ideal nunca iria transbordar e, para 
isto, uma política deveria ser colocada em prática para 
descarregar certa quantidade de água a cada ano. No 
entanto, se o fluxo de entrada do rio fosse muito baixo, 
o nível do reservatório se tornaria perigosamente 
baixo. O problema era: que política de descarte de água 
poderia ser definida, de modo a que o reservatório 
nunca transbordasse ou esvaziasse? 


Hurst estudou como o alcance do nível do reservatório oscilava em torno de seu nível 
médio. Se fluxos de entrada sucessivos eram aleatórias (ou seja, estatisticamente 
independentes) este intervalo - como acontece com o desvio-padrão no modelo de 
precificação de opções de Black-Scholes - iria aumentar ao longo do tempo de acordo com a 
sua raiz quadrada. Em seus estudos empíricos Hurst verificou que isto não acontecia surgindo 
então o conceito de memória de longo prazo (PETERS, 1991). 

Uma série de observações apresenta memória de longo prazo quando os valores 
observados em lags distantes são correlacionados entre si, ou ainda, se o efeito de um evento 
ocorrido em um instante pode ser detectado muitos lags depois. Teoricamente, Hurst (1951) 
foi o primeiro a observar a presença da memória longa ao constatar a presença de correlação 


entre um evento presente e eventos futuros, ou seja, uma oscilação causada hoje teria um 
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impacto longo. Posteriormente, Mandelbrot (1968) aprofundou o estudo sobre o tema de 
memória longa replicando seu conceito em demais áreas de estudos. 

Giraitis et. al. (2003) definem memória de longo prazo, no contexto de séries 
estacionárias em covariância, como aquelas que apresentam a propriedade de a função de 
autocovariância não ser absolutamente somável. A memória de longo prazo de uma variável 
X t é tipicamente caracterizada pelo decaimento hiperbólico da função autocovariância, 
definida por yu = Cov (Xk; Xo). Esse decaimento tem taxa ck 2d ~ ] , com 0 < d < 0,5. Caso X, 
tenha memória de longo prazo, tem-se, conforme Morettin (2011): 


lim 

M —>00 . 

j=-n 

A detecção de memória de longo prazo em uma série de observações indica, conforme 
Cajueiro et al. (2006), que existe um processo dinâmico envolvendo a variável observada e 
outros fatores, um caso oposto do movimento browniano. 

No caso da memória de curto prazo, a função autocovariância decai de forma 
exponencial, tendendo a zero rapidamente. Havendo apenas memória de curto prazo, a soma 
das autocovariâncias tem valor finito, um caso contrário ao da equação (1). 

3.3 O EXPOENTE DE HURST 

A correlação de séries temporais pode manifestar-se através de uma propriedade 
conhecida por persistência. É possível estudar esta propriedade em uma determinada série 
temporal, sensível às suas condições iniciais, com o auxílio do parâmetro expoente de Hurst, 
ou expoente H. 

De acordo com Feder (1988), as séries que possuem H > 0,5 são persistentes, ou seja, 


i =G ° 


(i) 


incrementos positivos no passado aumentam a probabilidade de incrementos positivos no 
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futuro. Edgar Palomas Molina (2002) classifica esta propriedade dizendo que se a tendência 
da série temporal tenha sido positiva no último período observado, é mais fácil que esta tenda 
a continuar sendo positiva do que negativa no período seguinte, e o contrário também segue o 
mesmo raciocínio. 

Já as séries que possuem H < 0,5 são antipersistentes, que mostram que para 
incrementos positivos no passado, incrementos negativos no futuro têm maior probabilidade. 
Fazendo uma adaptação do raciocínio de Edgar Molina acima pode-se dizer que se a 
tendência da série temporal tenha sido positiva no último período observado, é mais fácil que 
esta tenda a ser negativa do que positiva no período seguinte, e o contrário novamente segue o 
mesmo raciocínio. 

No caso do movimento browniano ordinário, o valor de H é 0,5. Isto implica na 
existência de um passeio aleatório (;randow walk). Neste caso constata-se a não existência de 
memória de longo prazo, já que, conforme Grimmentt e Stirzaker (1982), no passeio aleatório 
a movimentação depende apenas do passo anterior e não dos demais anteriores. Portanto, se 
não existe correlação entre eventos passados e futuros, rejeitando-se a hipótese de existir um 
processo de memória longa. 

O quadro seguinte sintetiza os principais pontos destes conceitos: 


Valores de H 

Evolução temporal 

Comportamento 

0,5 

Movimento browniano 

Passeio aleatório 

Entre 0 e 0,5 

Movimento browniano fracionário 

Memória Longa Anti Persistente 

Entre 0,5 e 1 

Movimento browniano fracionário 

Memória Longa Persistente 


Fonte: Peters (1991), elaboração própria 


Quadro 1 - Expoente de Hurst vs evolução temporal e comportamento. 

Durante seus estudos, para os períodos de cheia do Rio Nilo, Hurst encontrou H=0,72. 
Como seus resultados foram diferentes do esperado (H = 0,5), Hurst chegou a uma conclusão 


revolucionária para época de que as séries temporais destes fenômenos eram auto- 
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correlacionadas. A diferença entre o valor teórico e empírico do parâmetro H ficou conhecida 
por fenômeno de Hurst (PETERS, 1991). 

As séries persistentes e antipersistentes possuem memória de longo prazo, enquanto 
que as não-auto-correlacionadas são sem memória (PETERS, 1991). Para melhor ilustrar a 
propriedade persistência, a Figura 1 mostra o comportamento de três perfis distintos 
associados à expoente de Hurst 0,1, 0,5 e 0,9. 



Tempo 

Fonte: Mandelbrot (2004) 

Figura 1 - Movimento browniano ordinário e fracionário. 

Um bom exemplo desta relação de persistência é o movimento browniano fracionário 
que podemos relacionar, a título de ilustração, às flutuações do mercado de bitcoins, mais 
precisamente ao índice BPI. Se as flutuações deste mercado são persistentes, então haverá 
uma probabilidade maior de acentuadas valorizações e desvalorizações do valor do bitcoin em 
longo prazo, o que geraria dispersões grandes deste valor em torno da média à medida que o 
tempo fosse passando tornando-o atrativo para especuladores interessados em obter altas 
somas de lucros com estas dispersões. 
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Por outro lado, quando há antipersistência, as flutuações do BPI são pequenas ao 
longo prazo, comparando-as com as flutuações persistentes, por causa de menores dispersões 
em torno da média. Isto geraria um processo mais lento de valorização ou desvalorização do 
bitcoin com o tempo, tornando este mercado pouco atrativo para especuladores por causa das 
baixas somas de lucro proporcionadas pelo mesmo. 

Para flutuações do BPI brownianas teríamos um mercado com oscilações de preços 
independentes entre si, representando o caminho do meio entre a persistência e a 
antipersistência. 

Desta forma é possível concluir que o expoente de Hurst é um parâmetro capaz de 
descrever o comportamento de sistemas ao mesmo tempo complexos e caóticos, compostos 
por um número muito alto de variáveis que os influenciam e extremamente sensíveis a estas, e 
com isso poder ter uma noção de como emerge um comportamento global destes. Ou seja, o 
expoente de Hurst é capaz de evidenciar algumas características do mecanismo de auto- 
regulação destes sistemas (PETERS, 1991). 

3.4 O CÁLCULO DO EXPOENTE DE HURST 

Segundo Peters (1991) o expoente de Hurst é uma estatística robusta, na medida em 
que necessita de poucas hipóteses sobre o sistema a ser estudado, que permite a classificação 
de séries temporais em processos aleatórios e não aleatórios mesmo diante de séries aleatórias 
não normalmente distribuídas. É amplamente utilizado na estimação da memória de longo 
prazo, possuindo diversas formas de cálculo, cada uma com suas vantagens e desvantagens. 

Cajueiro et al. (2006) afirmam existir diversas formas de se calcular o expoente H. 
Entre as metodologias hoje são mais utilizadas a análise R/S (Rescaled Range) clássica, 
desenvolvida por Hurst (1951) e Mandelbrot (1972), a análise R/S modificada, de Lo (1991), 
a estimação do parâmetro de integração fracionária por regressão espectral, ou log- 
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periodograma, de Gewek e Porter-Hudak (1983), o estimador semi-paramétrico por log- 
periodograma, de Robinson (1995) e a análise DFA ( Detrended Fluctuation Analysis) 
desenvolvida independentemente nos trabalhos de Moreira, Silva, e Kamphorst (1994) e de 
Peng et al. (1994). Essas metodologias podem ser utilizadas sem que se tenha conhecimento 
dos fatores que atuam no processo gerador de preços, levando em conta apenas as séries de 
retornos ou volatilidades (ou resíduos dessas variáveis em relação a regressores) para as quais 
se deseja estimar a memória de longo prazo ou expoente H (CAJUEIRO et ah, 2006). 

Todas as metodologias descritas acima podem resultar em parâmetros de memória de 
longo prazo, seja o expoente H, ou o parâmetro de integração fracionária d , que podem ser 
utilizados na modelagem econométrica de processos com memória de longo prazo na média, 
como o ARFIMA (k, d, 1) ( auto regressive fractionally integrated moving average ), proposto 
por Granger e Joyeux (1980) e Hosking (1981). A relação entre a ordem de integração d do 
processo e o expoente H de Hurst é dada por (MORETTIN, 2011): 

d -H - 0,5 . ( 2 ) 

Neste trabalho serão utilizados os métodos acima citados para o cálculo do expoente H 
visando confirmar a HMF para o mercado de bitcoins. Porém o foco será na estatística R/S, 
pelo fato de ser a metodologia mais consistente, segundo Mandelbrot (2004), dentre as 
citadas. Conforme Cajueiro et al. (2006), esta técnica tem sido avaliada por diversos 
estudiosos, que buscam verificar a confiabilidade dessa medida. 

Dentre os pontos fortes, Cajueiro et al. (2006) destacam que, ao contrário da análise 
espectral, a análise R/S pode detectar ciclos não periódicos e ciclos com período igual ou 
maior que o período da amostra. Entre os pontos fracos, o método também detecta memória 
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curta, sem diferenciá-la da memória de longo prazo. Lo (1991) propôs uma estatística R/S 
modificada para solucionar esse problema. 

Por outro lado, Teverovsky, Taqqu e Willinger (1999) criticam a estatística R/S 
modificada de Lo e recomendam que esse método de análise seja empregado juntamente com 
outras ferramentas para a verificação desse tipo de previsibilidade. Além disso, conforme 
Cajueiro e Tabak (2007), a estatística modificada de Lo (1991) apresenta problemas de 
correlações de memória curta que interferem no cálculo de memória longa. 

A estatística R/S nada mais é do que a divisão da amplitude dos valores de uma série 
por seus respectivos desvios padrão de uma amostra. Entretanto, Hurst (1951) afirma que não 
se deve efetuar este cálculo apenas uma vez para uma série, ou seja, não basta dividir a 
amplitude total da série por seu desvio padrão. É preciso dividir a série em séries menores e, 
para cada uma destas subdivisões, calcular a amplitude e dividi-la pelos respectivos desvios 
padrão. Dessa forma será obtido um valor da estatística R/S para cada subdivisão da série. O 
autor observou uma relação entre a estatística R/S e o tamanho da amostra, de maneira que a 
estatística R/S é igual à metade do número de observações elevado ao expoente H. 

O método R/S é um método com bases empíricas. Surgiu com a tentativa de se 
resolver o problema de um reservatório finito submetido a um fluxo de entrada aleatório. A 
questão baseia-se na determinação do volume de um reservatório, conhecendo-se seu fluxo de 
entrada de água Ç(t), e considerando que seu fluxo de saída seja igual a média de ç (t), de 
maneira que este reservatório nunca seque ou transborde (FEDER, 1988). A Figura 2 ilustra 


este problema. 
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Fonte: FEDER, 1988 

Figura 2 - Esboço de um reservatório de água - fluxo de entrada ç(t) e um de saída <ç> T . 

Na Figura 2, a diferença acumulada entre o volume de água que entra e sai deste 
reservatório é X(t). O volume do reservatório para que ele não seque e nem transborde, R, é a 
diferença entre o máximo e o mínimo conteúdo do reservatório (FEDER, 1988, p. 151). 

Considerando um intervalo de tempo r , a média do fluxo de entrada será: 

• < 3 ) 

T t= i 

Seja X(t) a diferença acumulada entre o fluxo de entrada e sua média: 

= },(l<t <t) . (4) 

11 =1 T 

O máximo e o mínimo de (4) representam o volume máximo e mínimo de água que 
passará pelo reservatório no período r . Desta forma, o volume total do reservatório para que 
ele nunca seque nem transborde, deverá ser a diferença entre o máximo e mínimo de X , 
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R(r) = max X (t, r) la < r - min X(í, r)^ r . (5) 

Obviamente, R(z) depende do fluxo Ç(t), que por sua vez depende do período r . 

A Figura 3 mostra as grandezas X(t), ç(t) e R para o lago africano Albert. Na figura, o 
volume do reservatório para que ele não seque e nem transborde é indicado por R (FEDER, 
1988). 



Fonte: FEDER, 1988 

Figura 3 - Volume anual das águas do lago africano Albert. 


Hurst investigou os períodos de cheia do Rio Nilo considerando os reservatórios 
finitos: o Lago Vitória e Lago Albert, os maiores afluentes do Rio Nilo, e percebeu que R 
dependia do período r como uma lei de potência (FEDER, 1988), 



( 6 ) 
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na qual H é o expoente de Hurst e S é o desvio padrão do fluxo de entrada de águas 
nestes reservatórios definido por: 


S = 


(\ 

u 


T 


X 1/ 

\/2 


Iifw-<4> 2 

í=i J 


(7) 


Calcula-se a média dos valores (R/S), e se associa ao tamanho do bloco r: 


(s/s), = i-£(fl/s), . ,8) 

n /= i 

Dados os pares (R/S) T e x obtidos, executa-se a regressão por mínimos quadrados 
ordinários obtendo-se finalmente o expoente H: 

log(R/S) r = log C + Hlog t . (9) 

O expoente de Hurst é determinado com intervalo de confiança de 95% (CAJUEIRO 
et ai., 2006), estimados pelo método dos mínimos quadrados ordinários. 

3.5 Calculando o expoente de Hurst ao longo do tempo 

O surgimento de novos campos de estudo foi favorecido com a observação da variação 
do expoente de Hurst ao longo do tempo. Muniandy, Lim e Murugan (2001), puseram em 
prática esta nova forma de análise. 

O cálculo do expoente de Hurst ao longo do tempo consiste em uma análise não 
apenas pontual (um valor de H somente calculado), mas sim em uma sucessão de valores de 


H calculados ao longo do tempo. 
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Cajueiro et al. (2004) afirmam que o ideal é que seja utilizada uma janela móvel na 
qual contenha, aproximadamente, 504 observações. O autor justifica tal fato, por ser este um 
número que se aproxima à quantidade de dias úteis de dois anos, tempo este que considera 
suficiente para os cálculos de memória longa e, assim, se minimizariam os problemas de 
correlações de memória curta que eventualmente poderiam existir. 

Para o mercado de bitcoins, como o índice BPI é avaliado a cada 60 segundos, 24 
horas por dia e 7 dias por semana pode-se utilizar uma janela móvel bem menor que 500 
observações. Além disto, como os dados estão disponíveis apenas a partir de julho de 2010, 
um valor razoável do tamanho da janela seria 90 observações ou 3 meses. Já para a análise em 
diferentes escalas, o expoente H será calculado, para cada dia da séries, em relação aos 
últimos 3, 6 e 12 meses, ou 90, 180 e 360 observações respectivamente. Conforme Cajueiro et 
al. (2004), uma janela inferior a 90 dias para o cálculo do expoente H via metodologia R/S 
pode viesar o cálculo e por isto este número foi escolhido. 

A principal vantagem da análise da variação do expoente de Hurst ao longo do tempo 
é o fato, de que, de maneira clara pode-se observar padrões de comportamento e tendências 
das séries dos valores do expoente H. 

3.6 O EXPOENTE DE HURST E A EFICIÊNCIA DOS MERCADOS 

Confirme verificado na seção 3.3, o expoente H possui informações sobre a 
persistência. Além disto, é possível extrair informações do ponto de vista econômico, isto é, 
sobre o grau de eficiência de um determinado mercado. 

Observando a eficiência de mercados de bolsa de valores, Cajueiro e Tabak (2004) 
verificaram que esta está associada com H =0,5. Segundo eles, a velocidade na transmissão de 
informação em mercados eficientes é muito grande e há uma grande sincronia entre os 
agentes destes mercados que transmitem estas informações. 
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Desta forma, eles acharam razoável considerar as maiores economias do mundo 
globalizado como as mais eficientes. Com esta informação foi possível observar a ascensão de 
mercados emergentes. Para testar esta ascensão, eles selecionaram um banco de dados de 
cotações diárias de índices de Bolsas de Valores de alguns países emergentes e as do Japão e 
Estados Unidos entre janeiro de 1992 e dezembro de 2002. Para calcular o expoente de Hurst 
destas flutuações eles utilizaram o método R/S e chegaram à conclusão de que as flutuações 
dos índices das bolsas de valores do Japão e dos Estados Unidos tendiam a 0,5 enquanto que 
os países emergentes se distanciavam mais deste valor, tanto para cima quanto para baixo. 

Portanto é possível concluir que o grau ideal de eficiência dos mercados está associado 
a H = 0,5. Pode-se utilizar este parâmetro para observar quanto os mercados de um país 
emergente está distante do grau ideal de eficiência. Possivelmente o tamanho desta distância 
está associado a inúmeros fatores, tais como: nível de distribuição de renda, de produção 
industrial, de emprego, de produção agrícola e o investimento em educação da população 


nestes países. 
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4 METODOLOGIA 

Nesta seção é apresentada a metodologia utilizada para a realização desta pesquisa, 
abordando a fonte de dados adquiridos, os passos para análise e tratamento dos mesmos. Na 
sequência serão descritos os passos para avaliar a hipótese do passeio aleatório. Em seguida 
serão descritos os métodos para calcular o expoente H em toda a série analisada com o intuito 
de confirmar a presença de memória longa rejeitando assim a HME. Finalmente será 
detalhada a análise do expoente H tanto ao longo o tempo quanto em três diferentes escalas 
para cada dia da amostra a fim de verificar a possibilidade de previsão de eventos turbulentos. 

4.1 Base de dados 

Os dados utilizados neste trabalho são as taxas diárias de fechamento do BPI (índice 
de preços do mercado de bitcoins ), expressas em dólares, divulgadas pelo Coindesk, no 
período de 18 de julho de 2010 a 30 de setembro de 2014. Os dados foram obtidos da API do 
site do Coindesk. A Tabela 1 apresenta um resumo da série analisada. 


Tabela 1 Série analisada - dados diários - cotações de fechamento BPI. 


Série 

Moeda 

Período de/até 

Observações 

BPI 

US$ 

18/07/2010 a 30/09/2014 

1536 


Fonte: Coindesk, 2014, elaboração própria. 


Vale ressaltar que a escolha do BPI como proxy da análise do expoente de Hurst 
aplicado ao mercado de bitcoins se justifica uma vez que, conforme o capítulo 2, seção 2.8, 
representa os índices das quatro principais bolsas de bitcoins e pela sua relevância neste 
mercado (COINDESK, 2013). 

As Figuras 4 e 5 apresentam os gráficos do índice BPI e de seu logaritmo, 
respectivamente, durante o período analisado. 
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BPI - US$ 



2010 2011 2012 2013 2014 

Fonte: COINDESK, 2013, elaboração própria 

Figura 4 - Cotações diárias do BPI, em dólares. 


LOG BPI - US$ 



2010 2011 2012 2013 2014 

Fonte: COINDESK, 2013, elaboração própria. 

Figura 5 - Cotações diárias do logaritmo do BPI. 









47 


4.2 Análise e tratamento dos dados 

Com o intuito de se verificar a presença de memória longa persistente nas séries de 
fechamento do índice BPI, utilizaram-se dados de retornos diários, adquiridos a partir da série 
original, através da expressão R t = log (X t ) - log (X t .i), sendo X, e X t -i os fechamentos dos 
índices nas datas t e t - 1. 

A caracterização da série de log-retornos obtida acima é apresentada na Tabela 2. O 
gráfico correspondente a esses dados é apresentado na Figura 6. Para a utilização dos 
procedimentos estatísticos utilizou-se o software EViews. 


Tabela 2 Caracterização da série de log-retorno do BPI. 


Medida / Estatística 

Valor 

Média 

0,55E-2 

Mediana 

0,19E-2 

Máximo 

0,42 

Mínimo 

-0,49 

Desvio-padrão 

0,07 

Assimetria 

-0,33 

Curtose 

12,35 

Jarque-Bera 

5625,09 

Probabilidade 

l,27E-08 

Observações 

1536 


Fonte: COINDESK, 2013, elaboração própria. 
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4.3 Avaliação da hipótese de passeio aleatório 

Antes de verificar a Hipótese de Mercados Eficientes (HME) da série temporal dos 
log-retornos do BPI tem-se como primeiro passo testar a hipótese de que a série possa ser 
descrita como um processo do tipo random walk. 

De acordo com Santos (2013), a identificação de um processo do tipo random walk ou 
movimento browniano em séries temporais se dá pela verificação de algumas características, 
tais como a linearidade da série temporal, a distribuição de probabilidade dos log-retornos e, 
principalmente, a sua estacionariedade. 

Para avaliação da linearidade serão utilizados os testes de razão de variância (ou 
variance ratio tests) e também teste BDS (Brock, Dechert e Scheinkman, 1987) capaz de 
detectar diversas situações em que as variáveis não são independentes e identicamente 


Log-retornos BPI 



2010 2011 2012 2013 2014 

Fonte: COENDESK, 2013, elaboração própria. 

Figura 6 - Série dos log-retornos do BPI. 
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distribuídos (i.i.d.) além do caos determinístico, tais como a ausência de estacionariedade e a 
não-linearidade (SANTOS, 2013). 

Por meio de testes de correlação serial, de raiz unitária e estatísticas Portmanteau, será 
avaliada a estacionariedade da série, que, conforme Santos (2013), são bem representados 
pelos testes de Ljung-Box, de Phillips e Perron (1988) e o teste KPSS (Kwiatkowski, Phillips, 
Schmidt e Shin, 1992). 

Por último, a normalidade dos dados será avaliada investigando o coeficiente de 
curtose, K, indicativo do grau de achatamento da distribuição (SANTOS, 2013). 

Todos os procedimentos estatísticos e testes de hipóteses foram realizados através do 
software EViews. 

4.4 A CONFIRMAÇÃO DA PRESENÇA DE MEMÓRIA LONGA 

Conforme mencionado na seção 3.3, o expoente de Hurst é capaz de evidenciar a 
presença de memória longa persistente em séries temporais, rejeitando desta forma a Hipotése 
de Mercados Eficientes (HME). 

A análise R/S (Rescaled Range) clássica desenvolvida por Hurst (1951) e Mandelbrot 
(1972), já evidenciada e detalhada na seção 3.4, é, segundo Peters (1991), uma das mais 
difundidas técnicas para a detecção da fractalidade nas séries temporais. Ela é o foco principal 
deste trabalho, mas como a técnica tem sido avaliada por diversos estudiosos, que buscam 
verificar a confiabilidade dessa medida, outras metodologias serão utilizadas para fins de 
comparação de resultados. 

Dentre os estudos utilizando e comparando os principais métodos vale ressaltar o 
trabalho de Teverovsky, Taqqu e Willinger (1995) bem como o de Weron (2002). Neste 
último foram testadas as análises R/S clássica e R/S modificada de Lo (1991) bem como a 
DFA (Detrended Fluctuation Analysis) desenvolvida independentemente nos trabalhos de 
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Moreira, Silva, e Kamphorst (1994) e de Peng et al. (1994). Neste trabalho também foi 
utilizada a estimação do parâmetro de integração fracionária por regressão espectral, ou log- 
periodograma, de Gewek e Porter-Hudak (1983) mais conhecida como GPH. Os resultados 
evidenciam que a técnica DFA é suficientemente superior (WERON, 2002). Por outro lado, o 
trabalho de Kristoufek (2010) analisa as técnicas R/S e DFA e critica os resultados de Weron 
(2002), defendendo a R/S. 

Diante de diversas técnicas, cada uma com seus pontos fortes e fracos, o presente 
estudo irá, primeiramente, calcular o expoente H para toda a série temporal dos log-retomos 
do BPI, ou seja, no período de 18/07/2010 a 30/09/2014 contemplando 1536 observações. 
Foram utilizadas as metodologias R/S clássica e R/S modificada, bem como a GPH, o 
estimador semi-paramétrico por log-periodograma, de Robinson (1995) e a técnica DFA. 

Todos os procedimentos estatísticos e testes de hipóteses foram realizados através do o 
software R com a utilização dos pacotes nonlinearT series, fracdiff, fractal e pracma. O 
software EViews também foi utilizado para o desenvolvimento de um programa para o cálculo 
da estatística R/S disponível no apêndice. 

4.5 Avaliando o expoente de Hurst ao longo do tempo 

Conforme descrito na seção 3.5, o cálculo do expoente de Hurst ao longo do tempo irá 
analisar uma sucessão de valores de H calculados ao longo de uma janela móvel de 3 meses 
ou 90 dias deslocando-se de 10 em 10 dias. Conforme Cajueiro et al. (2004), uma janela 
inferior a 90 dias para o cálculo do expoente H via metodologia R/S pode viesar o cálculo e 
por isto este número foi escolhido. 

Para cada janela de 90 dias, calculou-se o expoente H através da técnica R/S clássica 
que, conforme Fo (1991), apesar dos pontos fortes deste método serem, até certo ponto, 
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passíveis de contestação, a estatística R/S clássica pode, reconhecidamente, identificar 
memória de longo prazo. 

Foram computados desta forma 145 valores de H ao longo de toda a série de 1536 dias 
através de um programa desenvolvido no software EViews disponível no apêndice. Se 
guardarmos o valor de cada expoente H calculado em um vetor de 145 posições temos, para 
cada posição K ( 0<K<146 ) deste vetor, o intervalo de dias considerados no cálculo dado por: 

10 * |X -1] +1 -> 10 * [£ -1] + 90 . (10) 

A fim de verificar a possibilidade de previsão de eventos turbulentos causadores das 
grandes distorções na série histórica dos log-retomos do BPI, analisou-se o efeito memória 
em três diferentes escalas de tempo para cada dia da amostra. De acordo com Mandelbrot 
(2004) e Feder (1988), a autossimilaridade indica que parte do objeto ou estrutura em análise 
é semelhante e está de alguma forma relacionada com o todo ao qual ela pertence. Ou seja, o 
expoente de Hurst calculado em diferentes escalas temporais deveria, teoricamente, ser 
semelhante para todas as escalas. 

Desta forma, uma segunda análise do expoente H ao longo do tempo é o seu cálculo 
em diferentes escalas, ou seja, avaliar seu valor para cada dia da série em relação aos últimos 
3, 6 e 12 meses, ou 90, 180 e 360 observações respectivamente. 

Novamente foi desenvolvido um terceiro programa no software EViews disponível no 
apêndice análogo ao anterior com particularidade de serem computados 3 vetores contendo os 
valores de H para cada dia da amostra, cada vetor representando cada uma das 3 escalas. 
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5 RESULTADOS 

5.1 Testes da hipótese do passeio aleatório 

Conforme metodologia descrita na seção 4.3 e seguindo um procedimento análogo ao 
realizado por Santos (2013), incialmente foi analisado o correlograma da série dos log- 
retomos do BPI e, avaliando as funções de autocorrelação e autocorrelação parcial, a série se 
demonstrou estacionária como pode ser verificado no anexo IV. 

Com o correlograma também foi possível realizar o teste de Ljung e Box modificado 
para verificação da ausência de autocorrelação nas observações componentes da série. Os 
resultados indicaram rejeição (exceto componentes 3 e 4) da hipótese nula de que não há 
autocorrelação, oferecendo, conforme Santos (2013) indícios de alguma relação de 
dependência, não sendo possível afirmar que a série dos log-retornos do BPI segue um 
processo aleatório com observações independentes. 

Em seguida foi realizado o teste de Phillips e Perron a fim de verificar a hipótese nula 
de que a série dos log-retornos do BPI apresentava raiz unitária, ou seja, 1(1). Os resultados 
indicaram a rejeição da hipótese nula, sendo necessário avaliar se a série seria integrada de 
ordem zero 1(0) ou fracionalmente integrada, ou seja, 0 < d < 1 (SANTOS, 2013). 

Desta forma foi realizado o teste KPSS para avaliar a hipótese nula que a série é 
estacionária apenas em torno de uma constante. Os resultados indicaram que não foi rejeitada 
a hipótese nula, havendo indícios de que a série seria estacionária. 

Mesmo rejeitando a hipótese nula da série ser 1(1) e não rejeitando a hipótese nula da 
série ser 1(0), sabe-se, conforme Santos (2013) da possibilidade de existência de memória 
longa e da baixa potência dos testes de raiz unitária diante de séries com esta característica. 

Realizou-se então o teste BDS, para verificar a hipótese de que as componentes da 
série dos log-retornos do BPI eram independentes e identicamente distribuídas (i.i.d.), o que 
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foi rejeitado, tanto com como sem bootsírap, indicando novamente a existência de alguma 
relação de dependência entre as observações. 

Além disto, avaliou-se a hipótese nula de que os dados são normalmente distribuídos 
por meio do teste de Jarque-Bera e da avaliação do coeficiente de assimetria e curtose da 
série. Avaliando os dados da tabela 2 da seção 4.2, da tabela 3 abaixo e pelo histograma nos 
anexos, nota-se que, novamente, a hipótese nula foi rejeitada, indicado um comportamento 
bem distante de uma distribuição gaussiana. 

Finalmente foi realizado o teste de razão de variância verificando hipótese nula de que 
as séries analisadas seguiam um passeio aleatório com erros heterocedásticos, variâncias 
viesadas e com bootstrap, pois o comportamento das probabilidades estava distante da 
normal. Os resultados evidenciaram comportamento distante da hipótese de passeio aleatório, 
uma vez que, conforme Santos (2013), era esperado uma razão de variância próxima de um. 

Todos os testes foram realizados através do software EViews e os resultados podem 
ser consultados nos anexos e foram resumidos da tabela 3 detalhada abaixo. Importante 
mencionar que o valor 0,00 para o p-valor é uma limitação do software devido ao número não 
suficiente de casas decimais. 


Tabela 3 Testes da hipótese do passeio aleatório da série de log-retornos do BPI. 


Teste 

Ho 

p-valor 

Resultado 

Ljung-Box 

Não há correlação 

0,00 

Ho rejeitadada 

Phillips-Perron 

Existe raiz unitária 

0,00 

Ho rejeitadada 

KPSS - intercepto 

Dados estacionários 

0,20 

Ho não rejeitadada 

BDS com bootstrap 

Dados i.i.d. 

0,00 

Ho rejeitadada 

Jarque-Bera 

Distribuição normal 

1.27E-08 

Ho rejeitadada 

Razão Variância com 

Dados martingales 

0,01 

Ho rejeitadada 

bootstrap 





Fonte: Resultados de teste pelo Eviews, elaboração própria. 
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5.2 Cálculo do expoente de Hurst pata toda a série 

Conforme resultados descritos na tabela 3, nota-se que existem fortes evidências de 
rejeição da hipótese do passeio aleatório para a série de log-retomos do BPI, fornecendo 
indícios de que a mesma poderia apresentar características de uma estrutura fractal, ou 
Hipótese dos Mercados Fractais (HMF). 

Desta forma, conforme mencionado na seção 4.4, a confirmação da presença de 
memória longa mediante a análise do expoente de Hurst foi fundamental para a identificação 
da HMF para toda a série de log-retomos do BPI no período de 18/07/2010 a 30/09/2014 
contemplando 1536 observações. 

Inicialmente foi utilizada a metodologia R/S clássica mediante programa desenvolvido 
no software Eviews disponível nos anexos. O valor encontrado para o expoente H foi de 0,65 
indicando comportamento de memória longa persistente. 

Em seguida, através do software R , foi calculado o valor do expoente H através das 
metodologias R/S clássica e R/S modificada, bem como a GPH, o estimador semi-paramétrico 
por log-periodograma, de Robinson (1995) e a técnica DFA. 

Todos os resultados podem ser consultados nos anexos e foram resumidos da tabela 4 
detalhada abaixo: 


Tabela 4 Valores estimados para o expoente Hurst. 


Metodologia 

H 

Erro Padrão 

R/S Clássica - Eviews 

0,65 

0,24E-2 

R/S Clássica - R 

0,66 

0,03 

R/S Modificada (Lo) 

0,66 

não disponível 

DFA 

0,64 

não disponível 

GPH 

0,53 

0,12 

Robinson 

0,61 

não disponível 


Fonte: Resultados da pesquisa via Eviews e R, elaboração própria. 
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Com relação aos intervalos de confiança, note-se que, conforme demonstrado nos 
anexos, as técnicas R/S Modificada de Lo, DFA e Robinson não computaram os valores do 
desvio padrão devido a uma limitação do pacote do Software R. De qualquer forma os 
resultados foram considerados satisfatórios confirmando o comportamento de memória longa 
persistente para a série analisada (H > 0,5) com exceção da metodologia GPH. 

Conforme mencionado na seção 4.4, cada técnica possui suas vantagens e 
desvantagens sendo importante a realização de testes comparando os resultados. Nota-se que 
a GPH possui valores de H abaixo de 0,5 no intervalo de confiança. O fato pelo qual a 
metodologia GPH ter rejeitado a HMF pode estar ligado, conforme Cajueiro et al. (2006), ao 
fato de as análises espectrais não conseguirem identificar a memória de longo prazo em séries 
temporais altamente afastadas da normalidade, com grande assimetria e curtose. 

5.3 CÁLCULO DO EXPOENTE DE HURST AO LONGO DO TEMPO 

De acordo com a metodologia apresentada na seção 4.5, o cálculo do expoente de 
Hurst ao longo do tempo foi realizado obtendo-se o gráfico da Figura 7 para os log-retomos 
do BPI, apresentando valores com intervalo de confiança de 95% para H, conforme o modelo 
explicitado pela equação (9). 

Nota-se na Figura 7 que o eixo das abscissas não contém datas, mas para cada posição 
K ( 0<K<146 ) do vetor dos 145 valores de H computados, o intervalo de dias considerados no 
cálculo é dado pela equação (10). Vale relembrar que foi utilizada uma janela móvel de 90 
dias deslocando-se de 10 em 10 dias. 

A série de expoentes de Hurst ao longo do tempo foi caracterizada estatisticamente e 
os resultados foram apresentados na Tabela 5. Foi calculada a estatística de Jarque-Bera eop- 
valor obtido (0,14) não permite rejeitar a hipótese nula de normalidade. 
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Série: LOG-Retornos do BPI (US$) 



índices do Vetor contendo os 145 Valores de H 
Fonte: Resultados da pesquisa via Eviews, elaboração própria. 

Figura 7 - Expoente de Hurst ao longo do tempo - log-retornos do BPI. 

A não rejeição da normalidade da série de expoentes H calculados indica que esses 
variam com o tempo de uma forma uniforme não podendo rejeitar a HMF (H > 0,5). 


Tabela 5 Caracterização da série de expoentes de Hurst calculados. 


Medida / Estatística 

Valor 

Média 

0,67 

Mediana 

0,68 

Máximo 

0,86 

Mínimo 

0,48 

Desvio-padrão 

0,08 

Assimetria 

-0,08 

Curtose 

2,21 

Jarque-Bera 

3,90 

Probabilidade 

0,14 

Observações 

145 


Fonte: Resultados da pesquisa via Eviews, elaboração própria. 
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A Figura 8 abaixo apresenta uma comparação das Figuras 6 e 7. Pode-se verificar que 
os ciclos de altos e baixos do expoente H se relacionam com a variação dos log-retomos: os 
menores valores de H coincidem com os períodos de pequenas variações nos log-retomos e 
vice-versa. Esta é uma observação importante, pois é um caso de arbitragem, ainda mais pelo 
fato da variação do H seguir a distribuição normal. Conforme Cajueiro et al. (2006), este é um 
dos pontos fortes da análise R/S que, ao contrário da análise espectral, pode detectar ciclos 
não periódicos e ciclos com período igual ou maior que o período da amostra. 


Log-retornos BPI 



2010 2011 2012 2013 2014 


Série: LOG-Retornos do BPI (US$) 



índices do Vetor contendo os 145 Valores de H 


Fonte: Resultados da pesquisa via Eviews, elaboração própria. 

Figura 8 - Comparação da série dos log-retornos do BPI vs expoente H ao longo do tempo. 
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Diante deste fato, o cálculo do expoente de Hurst em diferentes escalas de tempo foi 
realizado a fim de verificar a possibilidade de previsão de eventos turbulentos ou alterações 
bruscas na série de log-retornos do BPI. 

Para isso, o valor de H foi calculado em diferentes escalas de tempo, ou seja, avaliar 
seu valor para cada dia da série em relação aos últimos 3, 6 e 12 meses, ou 90, 180 e 360 
observações respectivamente, obtendo-se assim os gráficos da Figura 9 a seguir. 

Comparando os gráficos da Figura 9 bem como os da Figura 8, verifica-se que existem 
indícios da possibilidade de se prever uma alteração brusca na série de log-retomos do BPI: 
um período de alta no expoente H na escala de 90 dias antecede uma alta na escala de 180 
dias que por suas ver antecede um alta do H na escala de 360 dias. 

Reforçando o que já foi observado na Figura 8, esta alta do H repercute em alterações 
bruscas positivas na série de log-retornos do BPI. O raciocínio análogo pode ser realizado em 
períodos de baixa comprovando desta forma um comportamento de memória longa 


persistente. 
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Série: LOG-Retornos do BPI (US$) 



Série: LOG-Retornos do BPI (US$) 



Série: LOG-Retornos do BPI (US$) 



Fonte: Resultados da pesquisa via Eviews, elaboração própria. 

Figura 9 - Expoentes de Hurst ao longo do tempo observando últimos 360,180 e 90 dias 

respectivamente. 
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Vale a pena ressaltar que, de acordo com Cajueiro et al. (2006), a avaliação da análise 
R/S clássica como método de estimação de memória de longo prazo envolve incertezas no 
cálculo do valor do expoente de Hurst. 

Uma primeira fonte de incerteza se relaciona ao fato do cálculo via técnica R/S 
clássica ser assintótico, o que resulta em erro, o qual será maior quanto menor for o tamanho 
da amostra submetida à análise. Em segundo lugar, a estatística R/S clássica não tem 
distribuição conhecida sob a hipótese nula de não haver memória de longo prazo, de acordo 
com Cajueiro et al. (2006). Por último, dado que o valor de H é obtido por meio de uma 
regressão, calcula-se também o intervalo de confiança de 95% para esse valor, ou seja, com 
uma incerteza. 
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6 CONCLUSÕES E CONSIDERAÇÕES FINAIS 

Nos mercados de capitais, a corrente teórica de eficiência de mercado pressupõe que 
os preços variam aleatoriamente, sendo impossível identificar padrões ou fazer previsões. Já 
existem muitos estudos evidenciando que Hipótese dos Mercados Eficientes (HME) não é 
válida e, por sua vez, pressupõem que existem padrões de comportamento nos preços, que 
podem ser identificados, tal como a existência do processo de memória longa persistente. 

Diante disto, este trabalho analisou o comportamento do principal índice do mercado 
de bitcoins (BPI) no período julho de 2010 a setembro de 2014. Para testar a eficiência deste 
mercado ou evidenciar a existência da memória longa, foram realizados diversos testes da 
hipótese do passeio aleatório bem como realizado o cálculo do expoente de Hurst via 
determinadas metodologias tanto em todo o período da série quanto ao longo do tempo. 

As principais conclusões obtidas foram a não constatação de comportamentos 
aleatórios e que o comportamento de memória longa persistente foi verificado e comprovado 
para a série de log-retornos do BPI, não somente em todo o período analisado mas também ao 
longo do tempo. Além disso, a análise em diferentes escalas indicou a possibilidade de 
previsão de eventos turbulentos no índice. Finalmente foi possível comprovar a Hipótese de 
Mercados Fractais (HMF) para a série histórica de retornos do BPI no período de julho de 
2010 a setembro de 2014. 

Como sugestão para trabalhos futuros, recomenda-se o estudo não somente dos log- 
retomos mas também das volatilidades possibilitando uma análise mais criteriosa do risco. 
Além disto, uma análise em tempo real do expoente de Hurst, em diferentes escalas, através 
do framework Ruby on Rails, cuja API já está disponível para o BPI, permitirá a análise em 
diversas escalas de tempo possibilitando a previsão de eventos turbulentos no índice com uma 
maior precisão. 
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ANEXO I - PROGRAMA EVIEWS R/S PARA CÁLCULO DO EXPOENTE DE HURST 

subroutine hurst(series y) 
series ym 
!n = @obs(y) 
series RS_sum 
series RS_ave 
for !i=2 to !n 
!count = 1 
series R 

for !j = 1 to !n-!i+l step !i 
series cumser 
smpl !j !i+!j-1 
!S = @stdev(y) 
series ym = y - @mean(y) 
cumser = @cumsum(ym) 

R(!count) = (@max(cumser) - @min(cumser))/!S 
!count = !count+l 
next 

smpl 1 !count-l 
RS_sum = @cumsum(R) 

RS_ave(!i) = RS_sum(!count-1)/(!count-1) 

next 

smpl @all 

!start = @floor(10 a (S loglO (!n)*.25)) 

!end = @floor(10 a (@ logl0 (!n)*.75)) 
smpl !start !end 

equation eql.ls QloglO (RS_ave) c QloglO (@trend(0)) 
vector(l) hurst = eql.@coefs(2) 
smpl @all 

delete cumser r rs_ave rs_sum ym 
endsub 
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ANEXO II - PROGRAMA EVIEWS R/S PARA CÁLCULO DO EXPOENTE DE HURST AO LONGO DO 
TEMPO 

subroutine hurst_var(series z) 
smpl @all 

delete hurst3 eq3 dlhurst3 cumser r rs_ave rs_sum ym y 
!m = @obs (z) 

!tam = 90 
!int = 10 
!count2 = 1 

for !k = 1 to (@floor((!m-!tam)/!int)*!int)+1 step !int 
smpl !k !tam+!k-l 
series y = z 
!n = @obs(y) 
series RS_sum 
series RS_ave 
for !i=2 to !n 
!count = 1 
series R 

for !j = 1 to !n-!i+l step !i 
series cumser 
smpl !j + !k—1 !i+!j-l+!k-l 

!S = @stdev(y) 
series ym = y - @mean(y) 
cumser = @cumsum(ym) 

R(!count) = (@max(cumser) - @min(cumser))/!S 
!count = !count+l 

next 

smpl 1 !count-l 
RS_sum = @cumsum(R) 

RS_ave(!i+!k-1) = RS_sum(!count-1)/(!count-1) 

next 

smpl @all 
!start = !k 
!end = !tam+!k-1 
smpl !start !end 
!x = !k-1 

equation eq3.1s SloglO (RS_ave) c SloglO ( @trend(!x) ) 
series hurst3 

hurst3(!count2) = eq3.@coefs(2) 

!count2 = !count2+l 

delete cumser r rs_ave rs_sum ym y 

next 

smpl 1 !count2-l 
plot (hurst3) 
smpl @all 
endsub 

call hurst var (dlindex) 
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ANEXO III - PROGRAMA EVIEWS R/S PARA CÁLCULO DO EXPOENTE DE HURST EM 3 
DIFERENTES ESCALAS AO LONGO DO TEMPO 

subroutine hurst_var_scale(series z) 
smpl @all 

delete hurst4 eq4 hurst5 eq5 hurstô eq6 cumser cumser2 cumser3 
r r2 r3 rs_ave rs_sum rs_ave2 rs_sum2 rs_ave3 rs_sum3 ym ym2 
ym3 y y2 y3 
!m = @obs (z) 

!stp = 1 
!intl = 90 
!int2 = 180 
!int3 = 360 
!count2 = 1 

for !k = !int3 to !m step !stp 

smpl !k-!int3+l !k 
series y = z 
!n = @obs(y) 
series RS_sum 
series RS_ave 
for !i=2 to !n 
!count = 1 
series R 

for !j = 1 to !n-!i+l step !i 
series cumser 

smpl !j+!k-!int3 !i+!j-1+!k-!int3 
!S = @stdev(y) 
series ym = y - @mean(y) 
cumser = @cumsum(ym) 

R(!count) = (@max(cumser) - @min(cumser))/!S 
!count = !count+l 

next 

smpl 1 !count-l 
RS_sum = @cumsum(R) 

RS_ave(!i+!k-!int3) = RS_sum(!count-1)/(!count-1) 

next 

smpl @all 

!start = !k-!int3+l 
!end = !k 
smpl !start !end 
!x = !k-!int3 

equation eq4.1s SloglO (RS_ave) c SloglO ( @trend(!x) ) 
series hurst4 

hurst4(!count2) = eq4.@coefs(2) 

smpl !k-!int2+l !k 
series y2 = z 
!n2 = @obs(y2) 
series RS sum2 
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series RS_ave2 
for !i=2 to !n2 
!count = 1 
series R2 

for !j = 1 to !n2-!i+l step !i 
series cumser2 

smpl !j + !k—!int2 !i+!j-1+!k-!int2 

!S2 = @stdev(y2) 
series ym2 = y2 - @mean(y2) 
cumser2 = @cumsum(ym2) 

R2(!count) = (@max(cumser2) - @min(cumser2))/!S2 
!count = !count+l 

next 

smpl 1 !count-l 
RS_sum2 = @cumsum(R2) 

RS_ave2(!i+!k-!int2) = RS_sum2(!count-1)/(!count-1) 

next 

smpl @all 

! start = !k-!int2+l 
!end = !k 
smpl !start !end 
!x2 = !k-!int2 

equation eq5.1s SloglO(RS_ave2) c SloglO( @trend(!x2) ) 
series hurst5 

hurst5 ( !count2) = eq5.@coefs(2) 

smpl !k-!intl+1 !k 
series y3 = z 
!n3 = @obs(y3) 
series RS_sum3 
series RS_ave3 
for !i=2 to !n3 
!count = 1 
series R3 

for !j = 1 to !n3-!i+l step !i 
series cumser3 

smpl !j+!k-!intl !i+!j-1+!k-!int1 
!S 3 = @stdev(y3) 
series ym3 = y3 - @mean(y3) 
cumser3 = @cumsum(ym3) 

R3(!count) = (@max(cumser3) - @min(cumser3))/!S3 
!count = !count+l 

next 

smpl 1 ! count-1 
RS_sum3 = @cumsum(R3) 

RS_ave3(!i+!k-!int1) = RS_sum3(!count-1)/(!count-1) 

next 

smpl @all 

!start = !k-!intl+1 

!end = !k 

smpl !start !end 



!x3 = !k-!intl 

equation eq6.1s SloglO (RS_ave3) c @loglO( @trend(!x3 
series hurst6 

hurst6(!count2) = eq6.@coefs(2) 

!count2 = !count2+l 

delete cumser cumser2 cumser3 r r2 r3 rs_ave rs_sum 
rs_ave2 rs_sum2 rs_ave3 rs_sum3 ym ym2 ym3 y y2 y3 
next 

smpl 1 !count2-l 
plot (hurst4) 
plot (hurst5) 
plot (hurst6) 
smpl @all 
endsub 

call hurst var scale (dlindex) 
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ANEXO IV - RESULTADOS DOS TESTES DA HIPÓTESE DO PASSEIO ALEATÓRIO PELO EVIEWS 


CORRELOGRAMA DOS LOG-RETORNOS DO BPI 


Correloyrama dos LOG-Retornos do BPI 


Date: 12/09/14 Time: 21:27 
Sample: 7/18/2010 9/30/2014 
Included obseivations: 1535 
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RESULTADOS DO TESTE DE PHILLIPS-PERRON PARA OS LOG-RETORNOS DO BPI 


Phillips Perron Unit Root Test - LOG-Retornos do BPI 


Null Hypothesis: DLINDEX has a unit root 
Exogenous: Constant 

Bandwidth: 9 (Newey-Westautomatic) using Bartlett kernel 



Adj. t-Stat 

Prob* 

Phillips-Perron test statistic 

-37.35006 

0.0000 

Test criticai values: 1% levei 

5% levei 

10% levei 

-3.434406 

-2.863219 

-2.567712 


*MacKinnon (1996) one-sided p-values. 

Residual variance (no correction) 

HAC corrected variance (Bartlett kernel) 


0.004849 

0.005054 


Phillips-Perron Test Equation 

DependentVariable: D(DLINDEX) 

Method: Least Squares 

Date: 12/09/14 Time: 21:31 

Sample (adjusted): 7/20/2010 9/30/2014 

Included observations: 1534 after adjustments 

Variable 

Coefficient 

Std. Error 

t-Statistic 

Prob. 

DLINDEX(-1) 

-0.951745 

0.025513 

-37.30471 

0.0000 

C 

0.005259 

0.001785 

2.947143 

0.0033 

R-squared 

0.475996 

Mean dependentvar 

5.89E-05 

Adjusted R-squared 

0.475654 

S.D. dependentvar 

0.096229 

S.E. of regression 

0.069681 

Akaike info criterion 

-2.488475 

Sum squared resid 

7.438536 

Schwarz criterion 

-2.481519 

Log likelihood 

1910.661 

Hannan-Quinn criter. 

-2.485887 

F-statistic 

1391.641 

Durbin-Watson stat 

1.997520 

Prob(F-statistic) 

0.000000 
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RESULTADOS DO TESTE KPSS PARA OS LOG-RETORNOS DO BPI 


KPSS Unít Root Test - LOG-Retornos do BPI 


Null Hypothesis: DLINDEX is stationary 
Exogenous: Constaní 

Bandwidth: 11 (Newey-West automatic) using Bartlett kernel 


LM-Stat. 

Kwiatkowski-Phillips-Schmidt-Shin test statistic 

0.199944 

Asymptotic criticai values*: 1 % levei 

5% levei 

10% levei 

0.739000 

0.463000 

0.347000 

*Kwiatkowski-Phillips-Schmidt-Shin (1992, Table 1) 

Residual variance (no correction) 

HAC corrected variance (Bartlett kernel) 

0.004860 

0.005749 


KPSS Test Equation 

DependentVariable: DLINDEX 

Method: LeastSquares 

Date: 12/09/14 Time: 21:35 

Sample (adjusted): 7/19/2010 9/30/2014 

Included observations: 1535 after adjustments 


Variable 

Coefficient 

Std. Error t-Statistic 

Prob. 

C 

0.005480 

0.001780 3.078906 

0.0021 

R-squared 

0.000000 

Mean dependentvar 

0.005480 

Adjusted R-squared 

0.000000 

S.D. dependentvar 

0.069737 

S.E. of regression 

0.069737 

Akaike info criterion 

-2.487523 

Sum squared resid 

7.460202 

Schwarz criterion 

-2.484047 

Log likelihood 

1910.174 

Hannan-Quinn criter. 

-2.486230 

Durbin-Watson stat 

1.902840 
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RESULTADOS DO TESTE BDS PARA OS LOG-RETORNOS DO BPI 


BDS Test - LOG-Retornos do BPI 


BDS Testfor DLINDEX 
Date: 12/09/14 Time: 21:41 
Sample: 7/18/2010 9/30/2014 
Included observations: 1536 


Dimension 

BDS Statistic 

Std. Error 

z-Statistic 

Normal Prob. Bootstrap Pr... 

2 

0.049785 

0.003132 

15.89578 

0.0000 

0.0000 

3 

0.095292 

0.004988 

19.10277 

0.0000 

0.0000 

4 

0.124044 

0.005957 

20.82292 

0.0000 

0.0000 

5 

0.141381 

0.006228 

22.69920 

0.0000 

0.0000 

6 

0.149242 

0.006027 

24.76351 

0.0000 

0.0000 

7 

0.150198 

0.005542 

27.10078 

0.0000 

0.0000 

8 

0.146321 

0.004916 

29.76185 

0.0000 

0.0000 

9 

0.140468 

0.004246 

33.08320 

0.0000 

0.0000 

10 

0.133411 

0.003592 

37.14468 

0.0000 

0.0000 

Raw epsilon 


0.074357 




Pairs within epsilon 

1657365. 

V-Statistic 

0.703398 


Triples within epsilon 

2.01 E+09 

V-Statistic 

0.556103 


Dimension 

C(m,n) 

c(m,n) 

C(1 ,n-(m-1)) c(1,n-(m-1)) 

c(1,n-(m-1)) A k 

2 

639707.0 

0.544056 

826647.0 

0.703044 

0.494271 

3 

519655.0 

0.442531 

825367.0 

0.702872 

0.347239 

4 

431800.0 

0.368196 

824364.0 

0.702935 

0.244152 

5 

366358.0 

0.312802 

823094.0 

0.702769 

0.171421 

6 

315291.0 

0.269552 

821837.0 

0.70261 4 

0.120310 

7 

274105.0 

0.234648 

820650.0 

0.702517 

0.084450 

8 

239828.0 

0.205574 

819447.0 

0.702406 

0.059253 

9 

212461.0 

0.182354 

818955.0 

0.702905 

0.041886 

10 

189405.0 

0.162779 

817675.0 

0.702727 

0.029367 


RESULTADOS DO TESTE JARQUE-BERA PARA OS LOG-RETORNOS DO BPI 



Series: DLINDEX 

Sample 7/1 8/2010 9/30/2014 

Observations 1535 

Mean 

0.005480 

Median 

0.001865 

Maximum 

0.424580 

Minimum 

-0.491440 

Std. Dev. 

0.069737 

Skewness 

-0.333840 

Kurtosis 

12.35433 

Jarque-Bera 

5625.089 

Probability 

0.000000 
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RESULTADOS DO TESTE RAZÃO DE VARIÂNCIA PARA OS LOG-RETORNOS DO BPI 


Variance Ratio Test - LOG Retornos do BPI 

Null Hypothesis: Cumulated DLINDEX is a martingale 
Date: 12/09/14 Time: 22:07 
Sample: 7/18/2010 9/30/2014 
Included observations: 1535 (after adjustments) 

Heteroskedasticity robust standard error estimates 
Use biased variance estimates 
User-specified lags: 2 4 8 16 32 64 128 256 512 
Test probabilities computed using wild bootstrap: dist=twopoint, 
reps=2000, rng=kn, seed=1477928523 


JointTests 

Value 

df 

Probabilitv 

Max Izl (atperiod 64) 

3.186470 

1535 

0.0095 

Individual Tests 




Period 

Var. Ratio 

Std. Error 

z-Statistic 

Probabilitv 

2 

1.047922 

0.050612 

0.946847 

0.3470 

4 

1.039621 

0.090538 

0.437622 

0.6805 

8 

1.102080 

0.134066 

0.761418 

0.4370 

16 

1.263380 

0.182798 

1.440828 

0.1445 

32 

1.629220 

0.239024 

2.632462 

0.0090 

64 

1.994302 

0.312039 

3.186470 

0.0040 

128 

1.944021 

0.405206 

2.329731 

0.0105 

256 

1.383434 

0.535925 

0.715462 

0.5840 

512 

0.470024 

0.709344 

-0.747135 

0.7145 


TestDetails (Mean= 0.00548030587094) 


Period 

Variance 

Var. Ratio 

Obs. 

1 

0.00486 

- 

1535 

2 

0.00509 

1.04792 

1534 

4 

0.00505 

1.03962 

1532 

8 

0.00536 

1.10208 

1528 

16 

0.00614 

1.26338 

1520 

32 

0 00792 

1.62922 

1504 

64 

0.00969 

1.99430 

1472 

128 

0.00945 

1.94402 

1408 

256 

0.00672 

1.38343 

1280 

512 

0.00228 

0.47002 

1024 


Variance Ratio Statistic for Cumulated DLINDEX with Robust ± 2*S.E. Bands 
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ANEXO V - RESULTADOS DO CÁLCULO DO EXPOENTE DE HURST PARA TODA A SÉRIE DE LOG- 
RETORNOS DO BPI 


CÁLCULO DO EXPOENTE DE HURST VIA PROGRAMA EVIEWS - R/S CLÁSSICA 


DependentVariable: @L0G1 0(RS_AVE) 
Method: LeastSquares 
Date: 12/09/14 Time: 23:25 
Sample: 6 1426 
Included observations: 1421 


Variable 

Coefficient 

Std. Error t-Statistic 

Prob. 

C 

-0.170376 

0.006608 -25.78293 

0.0000 

@LOG1 0(@TREND(0)) 

0.652366 

0.002395 272.4299 

0.0000 

R-squared 

0.981239 

Mean dependentvar 

1.610142 

Adjusted R-squared 

0.981226 

S.D. dependentvar 

0.268387 

S.E. of regression 

0.036774 

Akaike info criterion 

-3.766658 

Sum squared resid 

1.918928 

Schwarz criterion 

-3.759256 

Log likelihood 

2678.210 

Hannan-Quinn criter. 

-3.763893 

F-statistic 

74218.03 

Durbin-Watson stat 

0.023895 

Prob(F-statistic) 

0.000000 












77 


CÁLCULO DO EXPOENTE DE HURT VIA SOFTWARE R 


A— METODOLOGIA GPH 

> fdGPH(DLINDEX, bandw.exp =0.5) 

$d 

[1] 0.02678344 -> H = d + 0,5 -> H = 0,52678355 


$sd. 
[ 1 ] 


$sd.reg 

[1] 0.1304275 


as 

0.1194042 


B— METODOLOGIA DFA 

> DFA(DLINDEX, detrend="polyl", sum.order=l, 

overlap=0.5, scale .max=trunc(length(DLINDEX)/1.5),scale.min=4,s 
cale.ratio=l.5, verbose=TRUE) 

Detrending model: x ~ 1 + t 
<environment: 0x038bfc98> 


Processing 

sca 

le 

4 

Processing 

sca 

le 

6 

Processing 

sca 

le 

9 

Processing 

sca 

le 

13 

Processing 

sca 

le 

20 

Processing 

sca 

le 

30 

Processing 

sca 

le 

45 

Processing 

sca 

le 

68 

Processing 

sca 

le 

102 

Processing 

sca 

le 

153 

Processing 

sca 

le 

230 

Processing 

sca 

le 

345 

Processing 

sca 

le 

518 

Processing 

sca 

le 

778 


Detrended fluctuation analysis for DLINDEX 


H estimate 

0.6361731 

Domain 

Time 

Statistic 

RMSE 

Length of series 

1535 

Block detrending model 

x ~ 1 + t 

Block overlap fraction 

0.5 

Scale ratio 

1.5 

Preprocessing 

lst order 


cumulative 


summation 
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Scale 4.000000 6.000000 9.00000 13.00000 20.00000 30.00000 
45.00000 68.00000 102.00000 153.00000 230.00000 345.0000 
518.0000 778.0000 

RMSE 0.056734 0.070742 0.08933 0.11035 0.14072 0.18343 

0.24256 0.32319 0.42939 0.59083 0.76941 1.0123 

1.3124 1.5226 


C-- METODOLOGIA R/S MODIFICADA (LO) 

> hurstexp (DLINDEX, d = 40, display = TRUE) 
Simple R/S Hurst estimation: 0.5777483 

-> Corrected R over S Hurst exponent: 0.663043 

Empirical Hurst exponent: 0.5879145 

Corrected empirical Hurst exponent: 0.5946988 
Theoretical Hurst exponent: 0.5300474 


D— METODOLOGIA R/S CLÁSSICA 

> rsFit (DLINDEX, leveis = 50, minnpts = 3, cut.off = 10 A c(0.7, 
2.5), 

+ doplot = FALSE, trace = FALSE, title = NULL, description 

= NULL) 

Title: 

Hurst Exponent from R/S Method 
Call: 

rsFit(x = DLINDEX, leveis = 50, minnpts = 3, cut.off = 

10 A c (0.7, 

2.5), doplot = FALSE, trace = FALSE, title = NULL, 
description = NULL) 

Method: 

R/S Method 

Hurst Exponent: 

H beta 

0.6577913 0.6577913 


Hurst Exponent 


X 0 


Estima te 
. 6577913 


Diagnostic: 

Std.Err t 


0.03072763 21 


-value 
40716 1 


Pr(>iti: 

47 68 94e-24 


Parameter Settings: 

n leveis minnpts cut.offl cut.off2 

1535 50 3 5 316 
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Description: 

Thu Oct 23 00:50:02 2014 by user: usuário 


E-- METODOLOGIA ESTIMADOR SEMI-PARAMÉTRICO POR LOG- 
PERIODOGRAMA, DE ROBINSON (1995) 

> hurstSpec(DLINDEX, method="smoothed", freq.max=0.25, 
dc=FALSE, n.block=NULL,weight=function(DLINDEX) 
rep (1,length(DLINDEX)), fit=lm, sdf.method="wosa") 

Hurst coefficient via regression of nonparametric sdf estimate 
for DLINDEX 


H estimate : 0.6121722 

Domain : Frequency 

Statistic : mean of SDF blocks 

Length of series : 1535 

Spectral Density Function estimation for 

Length of series 

1535 

Sampling interval 

1 

Frequency resolution (Hz) 

0.0006514658 

Centered 

TRUE 

Recentered 

FALSE 

Single-sided 

TRUE 

Method 

Wosa 

Block size 

383 

Overlap 

0.5 

Number of blocks 

Taper: hanning 

Number of points: 383 
Number of tapers: 1 

7 


Normalized: TRUE 
















